THE EXCITATION SPECTRUM FOR 
WEAKLY INTERACTING BOSONS 



ROBERT SEIRINGER 



Abstract. We investigate the low energy excitation spectrum of a Bose gas 
with weak, long range repulsive interactions. In particular, we prove that the 
Bogoliubov spectrum of elementary excitations with linear dispersion relation 
for small momentum becomes exact in the mean-field limit. 



1. Introduction and Main Results 

Bogoliubov's seminal 1947 paper [1^ contains several important results concern- 
ing the low energy behavior of bosonic systems. Among its striking predictions is 
the fact that the excitation spectrum is made up of elementary excitations whose 
energy is linear in the momentum for small momentum. Bogoliubov's method is 
based on various approximations and crucially uses a formalism on Fock space that 
does not conserve particle number. Mathematically, the validity of his method has 
so far only been established for the ground state energy of certain systems, see 
[El [El [HI HI m [E]- In particular, there are no rigorous results on the low energy 
excitation spectrum of interacting Bose gases, with the notable exception of exactly 
solvable models in one dimension [71 (H [ini HI [15] . 

In this article, we shall prove the validity of Bogoliubov's approximation scheme 
for a Bose gas in arbitrary dimension in the mean-field (Hartree) limit, where 
the interaction strength is proportional to the inverse particle number, and its 
range extends over the whole system. In particular, we verify that the low energy 
excitation spectrum for such a system equals the sum of elementary excitations, 
as predicted by Bogoliubov. As a corollary, we observe that the lowest energy in 
the sector of total momentum P depends linearly on |P|, a property that is crucial 
for the superfluid behavior of the system. The mean-field limit has served as a 
convenient and instructive toy model for several aspects of bosonic systems over 
the years. We refer to [51 [5] for a review and further references. 

We consider a homogeneous system of TV > 2 bosons on the flat unit torus T'^, 
d > I. The bosons interact with a weak two-body interaction which we write 
for convenience as {N — 1)^^11(2:). We assume that v is positive, bounded, periodic 
(with period one), and v{—x) — v{x). We also assume that v is of positive type, i.e., 
it has only non-negative Fourier coefficients. With A denoting the usual Laplacian 
on T*^, the Hamiltonian equals 

i?Ar = - ^ Ai + j^—j v{x^ - Xj) 
i—1 i<j 
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in suitable units. It acts on L'^y^{T^^), the permutation-symmetric square in- 
tegrable functions of N variables Xi £ T"^. Let Eo{N) denote the ground state 
energy of Hn- The Bogoliubov approximation [UIIl] predicts that Eo{N) is close 
to ^Nd{0) + E^°s, where 

= E (l^'l' + ^(P^ ~ VWT2\Wm) ■ (1) 

The sum runs over p e {2ttZ)'^, and 

v{p) = / u(a;)e~*P^da; 

are the Fourier coefficients of v. Note that the sum above converges, since the 
summands behave like v{p)'^ /\p\'^ for large p. 

More importantly, the Bogoliubov approximation predicts that the excitation 
spectrum of Hn is made up of elementary excitations of momentum p with corre- 
sponding energy 

ep = y/\p\^ + 2\p\^dip). (2) 
One noteworthy feature of ([2]) is that it is linear in p for small p, in contrast to the 
case when interactions are absent. 

Our main results can be summarized as follows. 

THEOREM 1. The ground state energy EaiN) of Hn equals 

Eo{N) = ^d{0) + E''°s + o{N-^/^) , (3) 

with E^°^ defined in {Ip. Moreover, the spectrum of Hn — Eq{N) below an energy 
^ is equal to finite sums of the form 

E e,n, + 0(e/'N-'/') , (4) 

pe(2TrZ)<i\{0} 

where Cp is given in (0) and Up £ {0, 1,2,...} for all p ^ 0. 

The error term 0(A^^^/^) in ^ refers to an expression that is bounded, in 
absolute value, by a constant times N~^/'^ for large iV, where the constant depends 
only on the interaction potential v] likewise for the error term 0(^^/^A^^^/^) in 
Q. The dependence on v is rather complicated but can be deduced from our 
proof, which gives explicit bounds. Keeping track of this dependence allows to 
draw conclusions about the spectrum for large N even if v depends on N . 

In the case of a fixed, TV-independent w. Theorem [T] implies that the Bogoliubov 
approximation becomes exact in the mean-field (Hartree) limit. As long as ^ <C 
N^/^, each individual excitation energy ^ is of the form X^p^p^-p with error o(l). 
Moreover, as long as ^ ^ iV, it is of the form X]epf^p(l + o(l)), i.e., the error is 
small relative to the magnitude of the excitation energy. In other words, in the 
mean field limit the whole excitation spectrum with energy ^ ^ is given in 
terms of sums of elementary excitations. The condition ^ <C A can be expected 
to be optimal, since only under this condition a large fraction of the particles are 
guaranteed to occupy the zero momentum mode, one of the key assumptions in the 
Bogoliubov approximation. For excitation energies of the order A" and larger, the 
spectrum will not be composed of elementary excitations anymore but has a more 
complicated structure. 
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Our proof shows that for each value of the {up} there exists exactly one eigen- 
value of the form ([4]). Moreover, the eigenfunction corresponding to an eigenvalue 
with given {up} has total momentum '^pPrip. Given this fact we readily deduce 
the following corollary from Theorem [TJ 

Corollary 1. Let Ep{N) denote the ground state energy of Hn in the sector of 
total momentum P. We have 

Ep{N) ~ Eo{N) = uim V + O (iFl^/^AT-iM . (5) 

In particular, Ep{N) - Eo{N) > \P\ miup y^2d(p) + \p\^ + 0(|P|3/2jV-i/2). 

The linear dependence oi Ep{N) on |P| is of crucial importance for the superfluid 
behavior of the Bose gas; see, e.g., the detailed discussion in [3l. It is a result of 
the interactions among the particles. The expression ([5|) differs markedly from the 
corresponding result for an ideal, non-interacting gas, especially if w(0) is large. 

Finally, we note that under the unitary transformation U — exp(— ig^^^j^ Xj), 
q G (27rZ)'', the Hamiltonian transforms as 

U'^HnU ^ Hn + N\q\^ ~ 2qP , 

where P = —i X^jLi denotes the total momentum operator. Hence our results 
apply equally also to the parts of the spectrum of with excitation energies close 
to N\q\'^, corresponding to collective excitations where the particles move uniformly 
with momentum q. 

The remainder of this paper is devoted to the proof of Theorem [T] The main 
strategy is to compare Hn with Bogoliubov's approximate Hamiltonian. The latter 
has to be suitably modified to take particle number conservation into account. 



2. Preliminaries 

We denote by P the projection onto the constant function in L'^{T'^), and Q = 
1 — P. The operator that counts the number of particles outside the zero momentum 
mode will be denoted by , i.e., 

N 

1=1 

We shall also use the symbol T for the kinetic energy T = ~ X^iLi ^i- 

Lemma [T] below gives simple upper and lower bounds on the ground state energy 
Eq[N) of i?Ar, as well as an upper bound on the expectation value of T in a low 
energy state. 

Lemma 1. The ground state energy of Hj^ satisfies the bounds 

N N 
> EoiN) - -d{0) > (^;(0) - v{0)) . (6) 

Moreover, in any N -particle state * with (*|iJAr|^') < yw(0) + p we have 

i2nf < (* \T\ < ^7^^ («(0) - ^(0)) +f,. (7) 
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Proof. The upper bound to the ground state energy follows from using a constant 
trial function. Since v > 0, we have 



pG(27rZ)<J\{0} 

This inequality can be rewritten as 



N 



2 

> 0. 



E vix.-^j)>^m~Y''^o) (8) 

l<i<j<N 

and thus 

iV N 
Hn > Y^iO) + T- («(0) - ?r(0)) . 

The rest follows easily. □ 

In the following, we shall also need a bound on the expectation value of higher 
powers of . More precisely, we shall use the following lemma. 

Lemma 2. Let 5* be an N -particle wave function in the spectral subspace of Hn 
corresponding to energy E < Eq(N) + fi. Then 

(27r)2 (5- \N>T\ < (^v{0) + 

N 

(m + 3w(0) + d{0)) (2^ + viO) - d{0)) . 



2(iV- 1) 

In particular, {^\N-^T\'i') is bounded above by an expression depending only on 
fi and V but not on N. 

Proof. Since 4* is permutation symmetric, 

(* \N>T\ *) iV (* \QiS\ *) + \N> {Hn - EoiN) - i^) | ^) 

where S — Eo{N) + — {N — 1}^^ ^i<j ^(^« ~ ^j)- Using Schwarz's inequality, 
the last term can be bounded as 



\N> {Hn - Eo{N) - ^^)\^^/) < L \{N>y\ ^) 
We split S into two parts, S = Sa + Sb, with 

S, = Eo{N) + ^-j^ v{x,-x,) 

2<i<j<N 

and 



1/2 



N 

I 

Sb = - 



1 

— ^w(a;i -Xj). 



N - ^ 

Note that Sa does not depend on xi. Using the positivity of v{p) as in dH), but 
with N replaced by A^ — 1, as well as the upper bound on Eq{N) in ([5]), we see that 

Sa < ^{n + d{0) + v{0)) . 

In particular, this implies that 

N (* \QiSa\ < ^ + w(0) + f (0)) |7V> I . 
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To bound the contribution of St, we use 

-{■^IQlSbl"^) = IQMxi - X2)\ ^) = \QlQ2v{xi - X2)\ *) 

+ {'i'\QlP2v{xi-X2)P2\-i') 
+ {'^\QlP2vixi~X2)Q2\'f) ■ 

The second term on the right side is positive. For the first and the third, we use 
Schwarz's inequaUty and |jw||oo = v{0) to conclude 

(^-igi^,!*) < v{0){^\QiQ2\^y^^ + i;(0)(^'|Qi|^') . 

Since 

A^'(«'|QiQ2|«') < (^'|(iv>)'|«') 

we have thus shown that 

|iV>r| «') < i (m + w(0) + 3t;(0)) \N> I 

+ {v{0) + ^^){^\{N>f\^f\ 

Using N> < {2'k)-^T this yields 

{^\N>T\^) < (-^-^^^^ + ifi + 3v{0) +v{0)) \N>\^) . 
The result then follows from Lemma [T] □ 

3. The Bogoliubov Hamiltonian 

The main strategy in the proof of Theorem [1] is to compare the Hamiltonian 
Hn with the Bogoliubov Hamiltonian. We will use a slightly modified version of 
it, which is, in particular, particle number conserving. 

Let Qp and denote the usual creation and annihilation operators on Fock space 
for a particle with momentum p G {2tt'E)'^, satisfying the canonical commutation 
relations [up, aj] — Spq. For p ^ 0, let 

bp = -^ii= (9) 

and define the Bogoliubov Hamiltonian 



(10) 



Note that conserves particle number. We are interested in H^°^ in the sector 

of exactly N particles. Note also that v{—p) — v{p) for all p E (27rZ)'^. 
Let 

Ap = \p\^ + d{p) 

and 

Bp = dip) . 
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A simple computation (compare with 12, Thm. 6.3]) shows that 
Ap (fet 5p + 6lp6_p) + Bp (^blblp + bpb-j)j 

([b\, + apb-p) ( bp + QpfeLp ) ( + apbp\ (b^p + apb^) 
\-al ^ \-al 



\ [Ap - ^Al ~ Bl) {[bp, bl] + [6_p, blp]) , 



where 



ap = [Ap - ^Al - Bl) ifBp>0, = if = 0. 



Note that < ap < +v{p)). In particular, supp_^Q ap < 1, and a(p) 

v{p)/\p\'^ for large p. 
Define further 

_ &p + ap^Lp 

Cp — I 

,^l-a2 

for p ^ 0. What the above calculation shows is that 

.Bog^_iv-^. [&P>^J] + [&-P.^- 

with Cp defined in The commutators equal 



[Vbt]^"°°°~"^"^<^. (12) 

4. Proof of Theorem [T] Lower Bound 

Recall that P denotes the projection onto the constant function in ^^(T''), and 
Q = 1 — P. Denote the two-particle multiplication operator v(xi — X2) hy v for 
short. Using translation invariance and the Schwarz inequality 

(P (g) Q + Q «) P)vQ ®Q + Q® Qv{P <^ Q + Q <^ P) 

> -e{P (E)Q + Q(g) P)v{P ®Q + Q®P)- e^^Q ® QvQ ® Q 
(which follows from positivity of v) we conclude that 

v>P® PvP (E)P + P(E) PvQ (E)Q + Q(E) QvP (g) P 

+ (1 - e)(P (g}Q + Q® P)v{P ®Q + Q®P)~ e'^Q ® QvQ ® Q (13) 

for any e > 0. The last term can be bounded by u(0)(5 ® Q- In second quantized 
language, this means that Hm is bounded from below by the restriction of 



Pt^O 



2e(iV- 1) 
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to the iV-particle sector. Here, we use again the definition ([9]) of bp. From now on, 
we shaU work with operators on Fock space, but it is always understood that we 
are only concerned with the sector of N particles. 
Next, we observe that 

t ^-1, + , 



Moreover, 



J2d{p)blbp<d{o)j^N> 

p/0 



and hence 



where the Bogoliubov Hamiltonian was defined in (ITU|) and 

1 ^ ^ iV>(iV>-l) f^^^^ ^ v{0)\ ^ 27V - 1 ^ > 



Ee = T+ ^ —^[d(0) + ^-^]+ev(0) N> . (14) 

' N-1 2{N-1) V ^-^TV-l ^ ' 

In particular, using PT|) and (IT^ . 

where E^°^ is the Bogoliubov energy defined in ((!]). 

The last term on the right side of (ITSt is positive and can be dropped for a 
lower bound on the ground state energy of Hn- For the choice e = 0{N~'^^^), the 
expected value of E^ in the ground state of is bounded above by 0{N~^^'^), 
as the bounds in Lemma [1] and [2] show. This proves the desired lower bound on 
EoiN). 

To obtain lower bounds on excited eigenvalues, it remains to investigate the 
positive last term in ([T5|) . We do this via a unitary transformation. Let U — , 
where 

x = Y.Pp{bibi^-bpb^p) 

with /3p > determined by 

tanh(2/?p) = ap . 

Note that X is anti-hermitian and hence U is unitary. If qq and aj were replaced 
by — 1, U would be the usual Bogoliubov transformation. Our modified U has 
the advantage of being particle number conserving, however. The price to pay for 
this modification is that U'^aqU can not be calculated anymore so easily. 
A second order Taylor expansion yields 

/■* 

-tx „ tx ,w „ 1 I / „\^~sXw w „ ^^„sx, 



Jo 

for any t > 0. We compute 

2 



8 



ROBERT SEIRINGER 



and 



=: 4/32a, + J, . (16) 
For any t > we thus have 



^a,e*-^ = ag + - — -Pqal^al + [t - s)e-''^ {Ap^aq + J,) e'^'^ds . 

— 1 Jo 

Iterating this identity leads to 

q2 

f/ta,t/ = cosh(2/3,)a, + s\nh.(2Pq)a-qj:^-^ + Kg dg + Kg 



with 



2,Pg 



In particular, we see that 

U^alagU = d\dg + KlKg + d\Kg + Kldg 

<{l + X)dldg + {l + \-^)KlKg (17) 

for any A > 0. We further have 

. _ t ( «o4 1^ a^g-paLp ^qp (ao)^Q§ ^ 



4. alan iV^ 



Using Schwarz's inequality, we can bound K^Kg as 

/fjA, < S25!|^ ^' .-.x,.j,,.x »inM2Wl-^)) ,. . 

To get an upper bound on JgJg, we write Jg as the sum of two terms, Jg^^ + 
where J^^^ is the second term on the right side in the first line of and Jg^"* = 
APlag{{N - l)-2ag(aS)2 - 1). Using jt J, < 2 J^^^^4^^ + 2jf jf ^ as well as 

we obtain the bound 

JlJq<C^Pq^-^^^^ ' ^i-64^(sup/3,^ + E/3^) • (20) 

To proceed, we need an upper bound on e-"^(iV> + ifeJ^ for < s < 1. For 
this purpose, let us compute 



[X,N>]^-2Y,(3q{blblg+bgb^g 
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We have 

- ^^^^^^ -l)iN + l- N>){N + 2-N>) 

+ iN> + 1){N> + 2){N - N>){N - 1 - N>)] 

^16(E/^p)]^(^> + ir- (21) 

With the aid of Schwarz's inequahty, we obtain 

[X, {N> + 1)2] = {N> + 1)[X,N>] + [X,N>]{N> + 1) 
> -7]{N> + lf -rr^[X,N>f 

for any 77 > 0. In particular, 

[X, (iV> + l)2]>-C2(iV> + l)2 

with 



We conclude that 

e-*^(7V> + l)2e*-^ < {N> + 1)2 + C2 f e-^^(7V> + ife'^'ds 

Jo 

for any t > 0. Iterating this bound gives 

In combination, ([19]), and ([22]) yield the bound 



K^K ^ (^" + l)V ,,c, (cosh(2/3,)-l)2 
By combining dTT]) with (HH]) and ([23l) . we obtain 



N>T 

lo.p]U — I sup ■ 



\ 9 



(cosh(2/3,) -1)2^ 



(24) 



Note that Iqp/S^ ^ w((7)2/|q|2 for large g, hence the last sum is finite. Applying this 
bound to (1151). we have thus shown that 



Hn + E,.^ > -v{0) + — 
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where 

N>T ( ep \ 

E,_^x =E, + — — sup / — ^ 

+ (1 + ^ ) iv-1 ^(^1.^^ I^^l J' 

with defined in ([Til) . 

The desired lower bound now follows easily from the min-max principle, and 
the fact that the spectrum of epaj^Up equals CpUp, with rip S {0, 1,2,...} 
for all p G (27rZ)''. In fact, for any function 4' in the spectral subspace of H]s[ 
corresponding to energy E < Eo{N) + ^, we have 

(*|i^e,A|*) < O ((e + iV-i) e + e'^"' (e"' + A^i)) 



according to Lemmas [T] and [2] The choice e = 0{^S^/N) = A then leads to the 
conclusion that the spectrum iJjv below an energy i?o(A^)+^ is bounded from below 
by the corresponding spectrum of 

^v{Q) + E^°^ +J2ep alap - O (^3/2^-i/2) . 

This completes the proof of the lower bound. 

5. Proof of Theorem [T] Upper Bound 

We proceed in essentially the same way as in the lower bound. In analogy to 
(UHl), we have 



V < P(g} PvP (E)P + P(E) PvQ (E)Q + Q(E) QvP (E) P 

+ (1 + e)(P (g}Q + Q® P)v{P ®Q + Q®P) + {1 + e-^)Q ® QvC 
for any e > 0. Together with 

blbp > alap ~ ^ 
this implies the upper bound 

Hn < Y^{0) + H^°^ + F, 

N> ^,,2N-1^ , i,iV>(iV> - l)w(O) , , 



where 



From a lower bound to the commutator ([T^ . namely 

2aJao - a^ap - aLpU^p _ _^ 



[bp, bi] + [b^p, bu . ^ "-^"-^ > 2 



TV- 1 - TV _ 1 ' 

we get 

^"°^<^"°<l-2(^)+E^.4^.- (26) 



p#0 
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Finally, to investigate the last term on the right side of this expression, we proceed 
as in (fT7 |) -([24 | . with the obvious modifications to get an upper bound instead of a 
lower bound. In replacement of (jl8p we use 



The result is 



dU >c^c 1 °S iV>(iV>+l) 



P \ P / \p^0 P ^ 



T 



sup 



iV-1 W{l-al) 



for any A > 0. 

Altogether, this shows that 

Hn < ^v{0) + £;^°s + ^^C/t ep alap^ U + F,^x , 

with F,,x given by the sum of in (gl]), ^N>{N - iy^\E^°s\ from and 
the last three terms in (j27p . To complete the upper bound, we need a bound on 
UF^ xU'^. For this purpose, we find it convenient to first bound F^ \ by 

Fe.^ < C, + A-i) ^^-^ + {e + N'^) 

for an appropriate constant C3 > 0. What remains to be shown is that 

f/(r+ < e'^''(T+ 1)2 (28) 
for some constant C4 > 0. Given we obtain 

+ Cse^^ (^(e-i + A-i) + {e + N-^) . (29) 

The spectrum of the operator on the right side of this inequality has exactly the 
desired form. Given an eigenvalue of X^p^^o '^\^p '^'^^^ value ^, we choose e = 
OiVUN) = A to obtain §v{0) + E^°^ + ^ + 0{i^/^N-^/^) for the right side of 
(P^ . This gives the desired upper bound. 

It remains to prove (|28p . This can be done in essentially the same way as in the 
proof of (1221) • In fact, 

and hence, similarly to (l2Tt . 

{N> + 
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Recall that ^ v{qY /\q\'^ for large which implies the finiteness of the sum since 
V is bounded by assumption (and hence, in particular, square integrable on T"^). 
By Schwarz, 

[X, [T+lf] <r^{T + lf +r]~^[X,Tf 
for any r/ > 0. In particular, 

[X, [T+lf] < CiiT+lf 
for some C4 > 0. We conclude that 

e^^{T +lfe-^^ < (T +lf + C4 [ e"'^ (T + l)^e""'^ds 

Jo 

for any t > 0. Iterating this bound yields (f28| . This completes the proof of the 
upper bound, and hence the proof of Theorem [1] 
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